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ON THE USE OF APPROXIMATE ANALYTICAL SOLUTIONS 


IN SOLVING OPTIMUM TRAJECTORY PROBLEMS 

By Donald J. Jezewski 
Manned Spacecraft Center 


SUMMARY 

A method is presented for obtaining numerical solutions of nonlinear, 
optimum trajectory problems by using an approximate analytical solution in a 
convergence process. Convergence to specified terminal conditions is based on 
an analytic solution representing' a reasonable approximation to the nonlinear 
system of equations. 

A numerical example of three-dimensional, time- optimal trajectory solutions 
in vacuum with a constant thrust is used to illustrate the technique. A gravi- 
tational force varying as the inverse square of the radius is considered. The 
approximate analytic solution is obtained from a system of quasi- optimum traJi- 
scendental equations. This analytic solution is used to obtain the initial and 
succeeding approximations to the initial values of the Lagrange multipliers. 
Generally, only a few iterations are required to obtain solutions which satisfy 
specified boundary conditions. 


INTRODUCTION 



The requirement for solutions of nonlinear, optimum trajectory problems 
has steadily increased in recent years. The nonlinear differential equations 
are difficult to solve because of the deficiency of a process which would con- 
verge rapidly to the specified boundary conditions. The first order perturba- 
tion techniques used in convergence processes normally work well in the vicinity 
of the solution where the linearity assun 5 )tion is valid, but they are usually 
unpredictable when large errors are present. 

The present paper deals with a technique which uses an approximate 
analytical solution in a convergence process for the corresponding nonlinear 
system of equations (refs. 1 and 2). The basis for convergence is that the 
solution curves of both the approximate and nonlinear equations are similar so 
that an identical change in the initial values of the Lagrange multipliers in 
each will produce approximately the same change in the terminal boundary con- 
ditions. These terminal boundary errors can be nullified by generating correc- 
tions to the initial values of Lagrange multipliers from the analytic solution. 
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SYMBOLS 

unit vectors in r, 0^ and cp directions, respectively 

azimuth angle measured from 1^, deg 

terminal "boundary vector 

initial values of Lagrange multipliers 

analytic solution 
general nonlinear functions 

function defined "by equation (4) 
altitude, ft 
specific impulse, sec 

logical control switch 
mass of vehicle, slugs 

radius to vehicle from center of attracting "body, ft 
thrust, lb 
time, sec 

components of total velocity in 1 , 1^, and 1^ directions, 
respectively, ft/sec 

total velocity, ft/sec 

initial weight of vehicle, lb 

inertial reference frame, ft 
state vector 
control vector 

function defined by equation (AlO) 
function defined by equation (A12) 
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Subscripts : 

0,1 


I 

R 


flight- path angle, deg 
latitude angle, deg 
Lagrange multipliers 

•z 2 

gravitational constant, ft^/sec 
longitude angle, deg 

thrust pitch angle measured from local horizontal plane, deg 
thrust yaw angle measured from 1^, deg 

initial and final values of state and control variables 
referring to a corrected solution 
referring to an integrated solution 
reference solution 


Superscripts ; 

vector quantity 
approximate value of variables 
T transpose matrix 


Operators : 

( ‘ ) time differentiation 


ANALYSIS 


Statement of Problem 

Consider a system of first-order differential equations of the form 
f(x,5^a,t) = 0 X (t^^ given 


to ^ t ^ ti 


where x and a axe respectively the state and control vectors. Certain 
terminal houndary conditions on the state vector are desired of the form 


B^[x(ti),ti] = 0 (2) 

The problem is to determine the control vector of (t) which forces the state 
from its given initial values boundary conditions defined by 

equation (2) . 

In optimization problems, the control vector is obtained from an optimality 
condition given by 


bE/ba = 0 


(3) 


where 


^ T 
H = A • 


'(x^x,a,t^ 




The superscript T refers to the transpose matrix, 
multipliers governed by the differential equation 


^ 


^ . 0 

ox 


The A's are Lagrange 


(5) 


This system of nonlinear equation (equations (l) to (5)) has a solution x (t) 
for any initial values of the Lagrange multipliers (^o) * initial control 

vector is obtained from equation (5) by the choice of the multipliers 7 \ ^t^^ . 
Thus, the problem is reduced to one of determining the initial vector 7 ^ (t^^ . 


Approximate Solution 

Consider an approximate analytic solution of the described system of 
nonlinear equations such that the solution evaluated at the terminal boundary 
conditions is 


/ 


where the circumflex over the variable indicates an approximate solution and 
the subscript R refers to a reference value. This solution can generally be 
derived by applying some simplifying assumptions to the nonlinear equations. 

A 

The initial values of the Lagrange multiplier from the analytic solution 7 ^ 

represent a first approximation to the Lagrange multipliers of the nonlinear 
system of equations. Ideally, they can be solved explicitly althou^ in 
practice the equations are usually transcendental, and iterative techniques must 
be applied. 


Iterative Technique 

Figure 1 is a logical flow chart of the iterative technique used in the 
convergence process. A switch k governs the various phases of the solution. 
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B I Resulting integrated vector 
T Lagrange multiplier 

Reference Lagrange multiplier 
Xc Corrected Lagrange multiplier 

Figure 1. - The use of an analytic solution in a convergence technique 

With the desired terminal boundary conditions and the switch k equal to 
0, the analytic solution is solved for the initial values of the Lagrange 
multipliers denoted as the reference solution. If the nonlinear system 

of equations is integrated to a terminal cutoff using the initial yalues x^ 
and the resiilting terminal vector ^ should resemble the desired terminal 

boundary conditions. The error in these conditions is indicative of the accuracy 
of the approximate solution. Replacing the desired boundary conditions by those 
obtained from the integrated solution = B^, the analytic solution is resolved 
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for the initial vector ?v. The switch k is set equal to +1. In the analytic 
solution, a change in the initial vector of 


a" - ^ (7) 

has produced a change in the terminal conditions of 


- b" (8) 

Since the analytic solution is an approximation to the nonlinear system of 
equations, it is assumed that the negative of this correction should produce an 
approximate change in the terminal boundary conditions of the integrated solu- 
tion. Therefore, the corrected initial value of the Lagrange multipliers for 
the nonlinear system of equations is 

a; = ^ (9) 

The switch k for the next and all future passes is set equal to -1. The 
procedure is repeated with one exception. The correction to the reference 

vector^ 7^ is obtained by subtracting the multiplier from the analytic solu- 
tion A^ and the previous Integrated solution 7 \ . This correction is given 
by 

✓s /s 

aA" = a"- A"^ (10) 

The solution is continued in this manner until the terminal boundary conditions 
approach the desired conditions within some arbitrary tolerance. Battin and 
Gibson (refs. 3 and 4) have used a similar technique in solving circumlunar 
trajectories using a matched conic as the analytic solution. 


DISCUSSION OF RESULTS 


The example problem that is used to illustrate the solution is a time- 
optimum descent from lunar orbit. The problem is typical of a number of 
examples which were solved. Refer to appendix A for a description of the 
numerical problem and the results of the derivation of the analytic solution. 

An IBM 704 o electronic data processing machine was used for the computa- 
tions. The solution time varied with the individual problem, depending mainly 
on the number of iterations required. For the longitude- free solution pre- 
sented, each Iteration required approximately 6 seconds. 
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' ' The iteration data and terminal boundary conditions are listed in table I 
for a time- optimum descent from a lunar orbit to an altitude of 5000 feet. The 
initial multipliers — those listed in the first iteration — vere obtained by 
evaluating equations (A13) to (A18) of appendix A from results obtained by the 
analytic, quasi-optimum solution. The equations given in appendix B were inte- 
grated forward using these initial conditions. The independent variable, time, 
was chosen as a stopping condition for the integration, since a value of this 
quantity is obtained from the analytic solution. For each subsequent iteration, 
the terminal time was corrected in the same manner as the initial, values of the 
Lagrange multipliers. Listed below the multipliers in table I are the terminal 
conditions on the boundary vector. A comparison of these values with the desired 
terminal boundary conditions shows good agreement, with the largest apparent 
error being in altitude. (This error is misleading, for the problem is actually 
solved in terms of the magnitude of the radius vector and not the altitude as 
the data inq)ly. The error in this variable is less than 1 percent for this 
iteration.) The solution can be observed to converge very rapidly; each itera- 
tion, which required only 5 or 6 seconds to compute, reduced the errors in the 
terminal boundary vector and the initial values of the Lagrange inultipliers by 
approximately 1 magnitude. A comparison of the multipliers in the first and 
last iteration indicates the excellent approximation of the analytic solution 
and the acceptability of the convergence technique. Total solution time for 
this problem was less than 30 seconds. 
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TABLE I.- TIME OPTIMUM DESCENT FROM LUNAR ORBIT 
(a) Iterations 



1 

2 

3 

4 


0.0 

0.0 

0.0 

0.0 

4 

Ag-lO .... 

-17.16130 

-22.23967 

-22.73874 

-22.75913 

A ^-10 . . • . 

5 

11.65346 

10.65851 

- .26852 

- .17464 



1.0 

1.0 

1.0 

1.0 

V-io^ .... 

5 

-1.219591 

-.880586 

-.945053 

-.943845 

Ag'lO^ .... 

-1.174105 

- 1.153440 

-1.173056 

- 1.173146 

ft/sec . . 

678.5133 

89.3276 

99.7368 

100.0205 

r ^, deg .... 

- 86.16649 

-7.18572 

.26271 

.00499 

deg . . . 

-108.6837 

-87.3814 

-86.9300 

-86.9981 

ft .... 

-71 477.5 

5411.5 

5062.5 

5000.4 

©2^^ deg .... 

-.366338 

- . 346616 

- .351411 

-.351574 

deg .... 

-28.02573 

-28.03312 

-28.03926 

-28.03928 

S 0 O • • • • 

326.0405 

327.6104 

327.4008 

527.5864 
1 


(b) Boxindary conditions 


T/Wq = 0.4172; I^^ = 514 secj 


Boundary 

conditions 

V, ft/sec 

r, deg 

deg 

h, ft 

0, deg 

cp, deg 

Initial 

Terminal 

5685 

100 

0 

0 

-83.5807 

-86.9994 

50 000 
5 000 

-1.5521 

-0.5516 

-17.95 

Free 
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CONCLUDING REMAEKS 


A technique has heen presented for obtaining solutions of optimum trajec- 
tory problems by using an approximate analytic solution in a convergence process. 
Convergence is based on the analytic solution representing a reasonable approxi- 
mation to the nonlinear equations. A numerical example of time-optimum, 
constant-thrust trajectory solution in three dimensions has been used to illus- 
trate the technique. The analytic equations were obtained from a quasi- optimum 
solution in which the initial values of Lagrange multipliers were an excellent 
approximation. Generally, only a few iterations were required to converge to 
the exact solution. 

Numerical results for a longitude-free solution are presented for a descent 
from lunar orbit. Solution time on an IBM 70^0 digital conqjuter was approxi- 
mately 50 seconds for this case. 


Manned Spacecraft Center 
National Aeronautics 
Houston, Texas, 


and Space Administration 
April 5 , 1965 
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APPENDIX A 


DERIVATION OF APPROXIMATE LAGRAHJE MULTIPLIERS 


The mathematical model employed is that of a mass particle with three 
degrees of freedom referred to a set of rotating coordinates. The axis system 
and the associated notation are illustrated in figure 2. The nonlinear differ- 
ential equations for this model and their associated Euler- Lagrange equations 
are given in appendix B. The vehicle is assumed to have a constant thrust and 
mass flow rate. No aerodynamic forces are considered. The angles x ^ 

are control variables and may be chosen arbitrarily. The final time t^ was 

selected to be minimized since, under the above assumptions, this will yield a 
minimum fuel consumption. 
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Figure 2. - Coordinate system and angie definition 


If, in the equations of motion, 
the nonlinear terms are replaced by 
power series in time and the calculus 
of variation is applied to this lin- 
earized set of equations, the following 
optimality conditions result (ref. 2) . 

tan X = (Al) 

sin Y = (A2) 


where Y has been assumed to be small. 
The multipliers are of the form 


\ = ^1 

(A3) 

CVI 

o 

II 

CVJ 

r< 

(a4) 


(A5) 

\ - qt 

(a6) 

Aj = Cj - C^t 

(A7) 


(A8) 


where to are constants. If the terminal boundary condition on the 
longitude is allowed to be free, the multiplier is 0. This result also 
facilitates integrating the linearized equations of motion. 
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(A9) 


Evaluating equations (Al) and (A2) at the initial and final time and 
maJcing use of equations (A3) to (a 8), the following relationships among the 
constants C will result: 


= "tan Xq 

= ("tan Xq - tan X^)/t^ 

Cg/Ci^ = sin Yq 

Cj/C^ = > (sin - sin Y^)/t^ 

Using these relationships, the Lagrange multipliers take the form 



^5 = tan Xq - P^t 
Xg = sin Yq - p^t 


(AlO) 

(All) 

(A12) 

(A13) 

(Al4) 

(A15) 

(A16) 

(A17) 

(A18) 


where the constant has been arbitrarily set equal to unity since the 

Euler- Lagrange equations are homogenous in the multipliers. 

From the solution of the analytic equations in reference 2, the initial 
and final values of the controls x an<3. Y along with the value of the final 
time tj^ are obtained. With these results, equations (A13) to (A18) are 

evaluated at the initial time to obtain a first approximation to the initial 
values of the Lagrange multipliers for the nonlinear system of equations . 

A result which is immediately apparent from equations (A15) to (Al8) is 
that a reasonably good approximation to the multipliers can be obtained by 
simply guessing five new parameters: Xq^ X^, Yq^ Y^, and t^. A knowledge 

of this set of variables for a particxxlar problem wovild certainly be more 
extensive than the set of multipliers. Except for cases in which the multi- 
pliers are extremely sensitive to the terminal conditions, as in a deceleration 
problem, this method should prove fruitful. 



APPENDIX B 


CONSTRAINT AND EULER- LAGRAMJE EQUATIONS 


The mathematical model is a mass particle with three degrees of freedom 
referred to a set of rotating coordinates. Hie axis system and the associated 
notation are illustrated in figure 2, In spherical coordinates, the equations 
of motion are: 


r - r(9^ + cp^ cos^ ©) ~ ^ " n/r^ (Bl) 

r 9 + 2r 0 + r cp sin 0 cos 0 = T/m cos x sin Y (B2) 

r 9 cos 9- 2r 9 cp sin 0 + 2r cp cos 0 = T/m cos x cos Y (B3) 


m = m^ - mt 


m = constant 


to ^ t ^ ti 


J 


(b4) 


where T/m is the acceleration due to the thrust and x and Y are arbitrarily 
chosen control variables. Let 


V = r 

u = r cp cos 9 


(B5) 

(b6) 


w 


r 9 


(BT) 


The constraint functions are 

f ^ = r cp cos 0 - u = 0 

fg = r - V = 0 

f_ = r0-w = O 
5 

f|^ = u - u (v - w tan 0)/r - T/m cos x cos Y = 0 

* / 2 2 \ / 2 
f^ = V + ^u + w j^r - T/m sin x + M*/r = 0 

fg = w - ^vw + u^ tan 9^ jr - T/m cos x sin Y = 0 


(b8) 

(B9) 

(BIO) 

(Bll) 
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(B12) 

(B15) 



The Euler-Lagrange equations corresponding to the six first-order 
differential equations described in equations (b 8) to (Blj) are 


• 

= - >v^(v - w tan e)/r 
^2 “ (^1 ^ ^3 ~ ^ 


(Bl4) 

- - QJ-A) + ^g('w + u^ tan 

(B15) 

• 

Aj = - ^uAj^ tan 0 + vA^^y^r - u(wA|j_ - uA^^ 

2 

sec 0 

2 

r 

(B16) 

A|^ = - A^ + [Ai^.(v - w tan e) - 2u^A^ - Ag 

tan 0 ^ J y r 

(BIT) 

• 

A^ = - Ag + (Ai^ u + Ag v^r 


(B18) 

Ag = - A^ - ^A|^ u tan 0 + 2A^ ^ 


(B19) 

The optimality conditions are 



Ac cos Y 
tan X = ^ . 

\ 


(B20) 

tan Y = Ag / Ai^ 


(B21) 
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